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ABSTRACT 

A nine-year record of attacks to Sitka spruce, Picea sitchensis (Bong.) Carr., by the 
Sitka spruce weevil (=white pine weevil), Pissodes strobi (Peck), was analyzed to 
determine the probability of attack on a tree based on the length of its terminal leader. 
Equations describing the relationship were developed. Tall trees with long leaders had 
higher rates of attack than short trees with short leaders. 
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RESUME 

Les donn£es compiles pendant neuf ans sur les d6gats causes par le char anon du pin 
blanc (Pissodes strobi Peck.) a des epinettes de Sitka (Picea sitchensis Bong.) ont 6t6 
analyses afin de determiner les probabilities qu’un arbre soit attaqu6 d’apres la longueur 
de sa pousse apicale. Des equations refletant cette relation ont 6\6 a abor6es. Les arbres 
de haute taille portant de longues pousses apicales 6taient plus fr6quemment attaqu6s 
que les arbres de faible hauteur a pousses apicales courtes. 


INTRODUCTION 

The Sitka spruce weevil (=white pine weevil), Pissodes strobi (Peck), is the most damaging 
pest of Sitka spruce, Picea sitchensis (Bong.) Carr., in coastal British Columbia, Washington 
and Oregon. In early spring, adult weevils crawl or fly to the terminal leader of the previous 
season and females lay eggs in niches excavated under the bark. The larvae kill the leader by 
mining and consuming the phloem. Following an attack, the lateral branches from the whorl 
immediately below the damaged leader develop negative geotropism and assume a vertical 
position. This process usually results in the formation of crooks and forks at the point of injury 
(Silver 1968, McMullen 1976, Alfaro 1989). Repeatedly attacked trees are stunted and 
overtopped by competing vegetation; a severely attacked plantation may be worthless (Alfaro 
1982). 

Modem pest management is greatly assisted by computer models that integrate pest 
biological and epidemiological factors into stand growth dynamics. Of particular importance 
are the factors that determine whether a particular stand or an individual tree within a stand is 
attacked. Several factors determine the rate at which Sitka spruce is attacked by the Sitka 
spruce weevil. Because of a climate unfavorable to insect development, stands in coastal 
locations and on the northern extremes of the distribution of Sitka spruce are less susceptible 
than stands on inland or southern locations (McMullen 1976, Heppner and Wood 1984). 
Similarly, a lower susceptibility of trees growing under shade, apparently caused by an 
unsuitable microclimate, has been reported for the eastern host of P. strobi, the eastern white 
pine (Pinus strobus L.) (Graham 1918; Wallace and Sulb'van 1985), as well as for Sitka spmce 
(McLean 1989). 

Plantations in the most susceptible areas are infested when trees are about 4 to 5 years old; 
attack intensity (number of trees attacked/year) increases rapidly thereafter, reaching a 
maximum when the plantation is 10 to 30 years old (Alfaro 1982; McMullen et al. 1987). The 
overall susceptibility of a stand decreases as trees reach heights above 12m (Harris et al . 1968; 
McMullen 1976, McMullen and Condrashoff 1973). This makes the Sitka spmce weevil a pest 
of primarily young forests. Overhulser et al. (1972) indicated that weevil oviposition and 
emergence is lower in trees that had previous attacks, relative to trees that had never been 
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attacked. Based on this finding, they hypothesized that older plantations have lower attack 
rates because they have a small number of unattacked trees and hence do not produce large 
numbers of weevils. 

Graham (1951) studied the attack records for a Sitka spruce plantation established in 
1930-1932 at Green Timbers Nursery, in Surrey, British Columbia. He concluded that the 
frequency of P. strobi attacks per tree could not be attributed to chance alone, but that some 
trees had higher rates of repeated attack than others. Graham did not speculate on the reasons 
for this non-randomness. In 1960, G.T. Silver, formerly with the Canadian Forestry Service in 
Victoria, British Columbia, established four plots to study the biology of this insect and 
possible means of control. Silver(1968) analyzed data collected in these plots between 1960 
and 1963 and reported that the tallest trees in a stand, with the longest leaders, had higher rates 
of attack than shorter trees with short leaders. Gara et al. (1971) and VanderSar and Borden 
(1977) corroborated Silver’s (1968) findings. McMullen et al. (1987) and McLean (1989) 
present figures that display the relationship between proportion or trees attacked and their 
leader length. However, these authors did not develop mathematical equations to quantify the 
relationship. The pictorial relationship in McMullen et al. (1987) is one of the several 
relationships used in a comprehensive population dynamics model these authors report in the 
same paper. 

McMullen et al. (1987) used data collected by Silver between 1960 and 1963. Silver 
continued the measurement of his plots until 1968. In this paper I analyze Silver’s full 9-year 
record and develop equations that describe the probability of a tree being attacked based on 
leader length and on weevil population level. 


MATERIALS AND METHODS 


Silver’s 1960s Study 

Four plots were established in the fall of 1959 in an area of natural Sitka spruce regeneration 
which originated after clearcut logging, near Nitinat Lake on Vancouver Island. The plots were 
rectangular, had a combined area of 1 ha, and initially included 692 trees which were marked 
with metal tags. One plot was left untreated and die others treated in certain years with 
insecticides. At the time of establishment, average tree age and height were 7 years and 1.3 m, 
respectively. In the early spring of every year, from 1960 until 1968, each tree was examined 
and tree height, the length of all leaders (including multiple leaders), and attack status (i.e., 
whether it was attacked or not) were recorded. Since the examinations were conducted before 
growth started, they represented tree condition at the end of the previous growing season. 

Data Analysis 

The data used in this paper were obtained from Silver’s check plot; they therefore represent 
uncontrolled damage levels. This plot occupied 0.32 ha and initially included 231 trees. 

To describe the stand and trees being attacked each year, the distribution of tree heights and 
leader lengths were tested for normality using the Kolomogorov test (Stephens 1983). This test 
was also applied to tree heights and leader lengths sorted by attack status. The data for every 
year were tested (Student’s t-test) for significant differences in height or leader length between 
attacked and unattacked trees. 

A logistic model (Hamilton 1974) was fitted to the binary attack data to relate the 
probability of a tree being attacked to its leader length. Leader lengths were converted to 
percentiles (Mendenhall 1975) of the leader length distribution for each year to allow 
comparisons between years, which were independent of the mean value. This is important 
because leader length increases with tree height and age up to a maximum which varies by site 
quality. Therefore, a particular length of leader may be considered long or short depending on 
plantadon height and age at the time. 

The probability of attack might depend not only on leader length but also on the level of the 
weevil population in a particular year. For this reason, in addition to the logistic model, 
separate linear models were calculated to describe the relationship between the percentage of 
attack in trees sorted by leader length percentile class in years of different attack severity. 
Severity classes were as follows: 
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LIGHT: the percentage of trees attacked in the year was 10% or less; 
MODERATE: the percentage of trees attacked in the year was between 10% and 
20 %; 

SEVERE: the percentage of trees attacked in the year was greater than 20%. 

The percentage of trees attacked refers to the percent of the stand trees attacked. Trees with 
multiple leaders were considered attacked if they had an attack in at least one of the leaders. 


RESULTS 

Tree height at the end of the 1959 and 1968 growing seasons averaged 1.3 and 4.6 m, 
respectively; the trees grew an average 3.3 m in the period. Leader length increased linearly 
with tree height until the trees were 3 to 4 m tall; it then reached a plateau at about 60 cm of 
growth per year (Figure 1). Fewer than 20% of the trees were attacked in the years from 1959 
to 1961. The percentage of trees attacked reached a maximum of 36% in 1964; thereafter, the 
attack rate oscillated around 30% per year (Table 1). 

The tests of normality indicated that, in every year, tree heights were distributed normally. 
Separate analysis of tree height distribution by attack status (attacked versus non-attacked) 
indicated no significant departure from normality (Kolomogorov D statistic test, P> 0.05). The 
distributions of the lengths of all leaders departed significantly from the normal distribution in 
four of the nine years studied. The difference, however, was only minor (probability of the D 
statistic was barely significant) and consisted of an excess frequency in the larger length 
classes. The distribution of leader lengths in leaders sorted by attack status were generally 
normal (Figure 2) with only two years in each class where a marginal departure from normality 
occurred, again due to an excess number of long leaders. The distributions of tree heights and 
leader lengths in attacked and non-attacked trees overlapped widely. In 1962, for example, the 
weevil attacked leaders as short as 25 cm but declined to attack many trees with 50- to 80-cm 
leaders. 

Piss odes strobi preferred the tallest trees with the longest leaders in all years (Table 1). The 
difference in height and leader length between attacked and non-attacked trees averaged 0.8 m 
mid 17.1 cm, respectively, over all years. 



FIGURE 1. Average leader growth of Sitka spruce tabulated by height at the start of the season (trees grouped by 
1-m height class) 
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TABLE 1. Percentage of Sitka spruce trees attacked by the Sitka spruce weevil and height and leader length of 
attacked and unattacked trees. 


Year® 

Attack b 

(%) 

Mean tree height (m) 
attacked not-attacked 

Mean leader length (cm) 
attacked not-attacked 

1959 c 

18 

2.0 

1.4 ** d 

37.0 

22.0 ** 

1960 

9 

2.3 

1.8 ** 

48.0 

33.0 ** 

1961 

14 

2.9 

1.9 ** 

56.7 

33.6 ** 

1962 

28 

3.1 

2.2 ** 

58.6 

37.0 ** 

1963 

33 

3.5 

2.3 ** 

62.5 

44.0 ** 

1964 

36 

3.7 

2.6 ** 

62.6 

55.1 ns 

1965 

35 

3.9 

3.1 ** 

64.0 

42.3 ** 

1966 

29 

4.3 

3.6 * 

58.7 

39.9 ** 

1967 

29 

4.6 

3.9 * 

61.6 

48.9 ** 


“Year of leader formation 

b Attacks occurred 1 year after leader formation 
c Average age in 1959 was 7 years. 

Significance of difference between mean tree heights or leader lengths were tested by the Student t-test, and are 
indicated by *=P< 0.05, **=P< 0.01, and ns= P> 0.05. 



LEADER LENGTH ( cm ) 

FIGURE 2. Frequency distribution of Sitka spruce leader lengths of trees attacked and not attacked by the Sitka 
spruce weevil, Pissodes strobi, in a representative year (1962), near Nitinat Lake, Vancouver Island. 

The logistic model: 

P = {1+EXP (A + BX LLENGTH)}- 1 (1) 

where P is the probability of a tree being attacked in a particular year, and LLENGTH is the 
length of the leader as percentile of the leader length distribution for that year, and a = 2.816, 
and b = -3.118, was highly significant (F=273, P<0.01). However, although this model 
provided a good estimate of the average probability of attack in any year (Figure 3), 
considerable variation remained unexplained (correlation coefficient = 0.33). This variation is 
due, in part, to the fact that the probability of attack on a tree depends not only on the length of 
the leader, but also on the population level of the weevil: in years of low population, a tree with 
a given leader length will have a lower probability of attack than in years with a high 
population. 
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FIGURE 3. Percentage of attack (dots) and probability of attack (solid line) by Pissodes strobi weevils in Sitka 
spruce leaders sorted by leader length class. Lengths were expresed as percentiles of the leader length distribution 
for each year. The width of the percentile length class was 10%. These data were collected over 9 years near 
Nitinat Lake, Vancouver Island. 


The linear models for each severity class were highly significant (Fig. 4). The models were 
refitted to eliminate non-significant intercepts, yielding the functions: 

Population level LIGHT: 

P= 0.002 x LCLASS 
r^O.68, F= 23.8 

Population level MODERATE 
P= -0.121 + 0.006 x LCLASS 
r 2 = 0.77, F= 60.9 

Population level SEVERE: 

P= 0.006 x LCLASS 
r2= 0.71, F= 145.2 

where P is the probability of a tree being attacked in a particular year; LCLASS is the midpoint 
of leader length class (Length is expressed as a percentile of the leader length distribution. 
Length classes have a width of 10% of the total length distribution, with mid-points at 5%, 
15%, 25%, etc.); and LIGHT, MODERATE and SEVERE are population levels as defined in 
the Methods section. 


DISCUSSION 

The attack rates increased from less than 20% to a maximum between 30% and 40% of the 
trees per year. This epidemiological pattern is very similar to that reported by Alfaro (1982) for 
a severely attacked plantation at Green Timbers, near Surrey, British Columbia and by 
McMullen et al . (1987) in plantations in the Klanawa River area of Vancouver Island. Alfaro 
(1982) indicated that attacked trees take 2 or 3 years to develop a new leader. During this 
period, the tree is generally not available for re-attack, unless the tree has developed multiple 
tops or the tree is re-attacked in the intemode below the previous year’s attack. In the Prince 
George area. Cozens (1987) found that in interior spruce (Picea glauca x engelmanni) up to 
20% of the attacks were re-attacks on trees attacked in the previous year. Therefore, attack 
rates of 30 to 40%, such as the ones reported here, are probably near the maximum attack rate 
that a weevil population may sustain in a stand. Higher attack rates would deplete the 
oviposition sites and would lead to a reduction in the weevil population. 
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FIGURE 4. Percentage of attack (dots) and probability of attack (solid line) by Pissodes strobi weevils in Sitka 
spruce leaders sorted by leader length class. A) When attack intensity is LIGHT i.e. the percentage of trees 
attacked in a stand is 10% or less; B) When attack intensity is MODERATE i.e. the percentage of trees attacked in 
a stand is greater than 10% but less than or equals to 20%; C) When attack intensity is SEVERE i.e. the percentage 
of trees attacked in a stand is greater than 20%. Lengths were expresed as percentiles of the leader length 
distribution for each year. The width of the percentile length class was 10%. These data were collected over 9 years 
near Nitinat Lake, Vancouver Island. Dotted line represents the 95% confidence interval. 
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The fact that the frequency distributions of attacked and non-attacked leaders overlap 
widely (Figure 2) suggests that, in addition to leader length and population level, other factors 
determine whether a tree is attacked or not. My own unpublished observations indicate that the 
spatial distribution of the weevil population and of the attacks in young plantations is clumped. 
Attack intensity in trees within population clumps would be higher than that between clumps, 
forcing the weevils to attack smaller leaders in the clump areas and leaving relatively longer 
leaders unattacked in the areas between clumps. 

The preference of P. strobi for the longest leaders in the stand is an adaptation of 
significance to the survival of a weevil population. It ensures that the leaders with the 
maximum food supply (more phloem in longer leaders) will be colonized; offspring produc¬ 
tion per leader is therefore optimized. 

The equations developed here could form an integral part of a pest management model such 
as the one presented by McMullen et al. (1987), to predict the damage caused by the Sitka 
spruce weevil to Sitka spruce in British Columbia. 
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